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Abstract 



Wc prove global existence in time of solutions to relaxed conservative cross diffusion systems 
governed by nonlinear operators of the form u% —$ dtUi — A(ai(u)uj) where the u%,i = 1,...,I 
represent / density- functions, u is a spatially regularized form of (u\, ...,itj) and the nonlinearities 
at are merely assumed to be continuous and bounded from below. Existence of global weak solutions 
is obtained in any space dimension. Solutions are proved to be regular and unique when the a.; are 
locally Lipschitz continuous. 

ON 

00 , 

5-) ' 1 Introduction 

(^ , Introduced by Shigesada et al. [28], cross diffusion models try to represent the effect of the interaction 

CN ' between species through motion, and not only as usual through reaction. These models have been 

studied by Levin [18], Levin and Segel [17], Okubo [26], Mimura and Murray [23], Mimura and 

Kawasaki [22], Mimura and Yamaguti [24], Andreianov et al. [1], Bendahmane and Langlais [2] and 

^ | many other authors: a survey by A. Jiingel may be found in [12] for applications to population 

H ■ dynamics. In those references, a general system is the following: 

d t ui - A[u 1 (d 1 + duu\ + d 12 v%)] = ri(u 1 ,u 2 ), 

d t u 2 - A.[u 2 (d 2 + d 2 iu I [ + d 2 2U^)] =r 2 (u 1 ,u 2 ), (1) 

d n [ui(di + dnUi + ^12^2)] = d n [u 2 (d 2 + ^21^1 + ^22^2)] = 0. 

For the system (1) with p = 1 and Lotka-Volterra-type reaction, there exists a wide literature, studying 
specific cases of the system where an additional structure keeps it parabolic or with cross diffusion 
pressure only on one of the species (see e.g. Wang [29] and the many references therein, especially 
in the introduction). To our knowledge, the most general result on global weak solutions might be 
found in Chen and Jiingel ["] where the entropy structure of the model is used. For existence of 
classical solutions the reader might consult [29, 19] by Wang and Li-Zhao for instance. In population 
dynamics, one of the most interesting features of cross diffusion is its effect on steady states: cross 
diffusion pressure might help the appearance of nonconstant steady states when the reaction structure 
does not drive to segregation (see Iida-Mimura-Ninomyia [11] for instance). However, in these cases, 



the pattern formation relies on the reaction term (for instance, the convergence to homogeneous steady 
states in absence of reaction is proved in [7]). 

In [3], [15], the first author and collaborators introduced a relaxation of conservative cross diffusion 
systems, replacing 

dtUi — A[ai(u)ui] = 0, on (0, +oo) X U, Q C M. N , bounded, 

u = (ui,. . . ,Ul), 

d n [di(u)ui] = on (0, +oo) x dQ, u(0, •) = u° given, 

where aj : [0, oo) 7 — > [0,oo), by the following relaxed model: 

dtUi — A[di(u)ui] = 0, on (0, +oo) x fl, 

u = (ui,...ui), , 2 > 

Ui — 5iAiii = Ui, on (0, +oo) x Q,, Si > 0, 

d n Ui = d n Ui = on (0, +oo) x Q,, u(0, •) = u° given. 

This model was introduced in order to investigate the effect of non classical cross diffusion pressure 
on the segregative behavior (and ctj(-) is often truly nonlinear). One of the purposes was to drive 
spatial segregation only through motion. Its effects on the stability of the homogeneous equilibria is 
investigated in [15, 3, 16]. This relaxed version is also relevant in some applications: it takes into ac- 
count that the intensity of the underlying Brownian motion depends on the density of the population 
measured with a space scale Si and not exactly at the exact location x. It takes therefore into account 
the fact that a species can react to the presence of another species in a neighborhood. 
Models with nonlocal diffusion coefficients can be seen also in [4] (where the self-diffusion coefficients 
depend on the total population). Nonlocal reaction terms can also be considered ([5, 8, 25] for in- 
stance), but the goal of our model is more to create patterns only through motion. 

A first well-posedness result for the relaxed model was derived in [15], [3] in dimension N = 1,2 
and with some restrictions on the structure of the nonlinearities a, (basically, the ai are C 2 and have 
at most a polynomial growth in u). In this paper, we prove existence of solutions for this system 
in any dimension and for general nonlinearities ai, which are only assumed to be continuous and 
bounded from below. Weak solutions are obtained in general and they are proved to be strong and 
unique as soon as the aj are locally Lipschitz continuous. Some L 2 -estimates are exploited in the 
spirit of [27] to prove existence of weak solutions. A main point is that u is uniformly bounded in 
any dimension for these weak solutions. Next, one has to deal with parabolic operators of the form 
Ui —> dtUi — A (ai(u)ui): they are not of divergence form, but they are uniformly parabolic since di(u) 
is then bounded from above and from below. Using the C Q -theory for the duals of these operators, 
namely Ui — > dtUi—ai(u)AUi, in the spirit of Krylov-Safonov [13], [10] (see also the book by Lieberman 
[20]), we prove that u is even Holder-continuous. This provides continuous coefficients ai(u) for the 
above operators, and then, L p -estimates classically follow for the solution. When the a% are locally 
Lipschitz continuous, even dtUi, A (ai(u)ui) are proved to be in L p so that the solution is strong: 
moreover, weak solutions are then proved to be unique. 

Let us fix the notations and state the main result. We assume that Q C R is a bounded subset 
with a C 2 -boundary. The exterior normal derivative operator on dQ is denoted by d n . For all T > 0, 
we denote Q T = (0,T) xfi,E T = (0,T) x dn. For a G (0, 1], we denote 

C a (Q T ) = {ve L°°(Q T ); \\v\\P < +oo }, 



I (") 



\l°°(q t ) +sup 



\v{t,x) -v(s,y)\ 

\t — s\ + \x — y\ 2 ]~- 



(t,x),{s,y) e Q T 



We will at least assume that 



Vi = 1, ...,/, a« : [0,oo) — > [0,oo) is continuous and : inf Oj(r) > d > 0. 

re[0,oo) 7 

And we are given <5j G (0,oo), Vi = 1, .../. 



(3) 



Theorem 1.1 Assume (3) and u° = (ui,...,u®) G L°°(fi, [0, 00)) . Then, there exists a nonnegative 
solution u = (ui, ...,u/) to i/ie following problem: 



( VTE (0,oo),Vi = l,...,/, VpG [l,oo), 

Uj G L p (Q T )-Ui G C Q (Q T ) n L p (0,T; W 2 ' p (fi)) /or so?ne a G (0, 1], 
^a^KGL^T;^ 2 *^)), 
Ui(t) - A[/ ai(u)ui] = u° in Q T , 
un — 5iAui = Ui in Qt 

d n f/o ai(u)uij = = d n Ui on S T . 



//" moreover 



(4) 



(5) 



Vi = 1, ..., J, a, : [0,oo) —7- [0,oo) is locally Lipschitz continuous 
then, Vi = 1, ...,/, VT > 0, Vp G [1, 00), 

Ui G L°°(Q t ),Vt G (0,^,^1^, A(ai(fi)ui) G L p ((r,T) x Q) 



and dfUi — A(ai(u)ui) = Q,d n (a,i(u)ui) = in a pointwise sense. Finally, under assumption (5), 
solutions of (4) are unique. 



The paper is organized as follows. 
Section 2 first assumes that the nonlinearities at are also bounded from above. We prove existence of a 
weak solution to the system (4) by a standard Leray-Schauder fixed-point argument. The underlying 
space is an adequate subspace of L 2 (Qt) and the required compactness follows essentially from Lemma 
2.3. 

Section 3 is devoted to the proof of the L°°-estimate on u. Then, the assumption of the bound from 
above on the a» may be dropped. 

Section 4 exploits this L°°-estimate to prove that the weak solution is actually rather regular, and 
existence as stated in Theorem 1.1 follows. The C Q -theory for nondivergential parabolic operators is 
used there. An alternative more elementary proof of the regularity is also given when monotonicity 
properties hold for the aj together with locally Lipschitz continuity. 

The uniqueness stated in Theorem 1.1 is proved in Section 5. It is based on solving an original dual 
problem, interesting for itself. 

A short Section 6 indicates without proof a complementary approach which provides a constructive and 
alternative way of proving existence of a solution and which may be used to compute it numerically. 



2 Global existence when ai is bounded 

In this section, we first prove existence of weak-solutions of (4) on a given interval [0, T] when, besides 
(3), the nonlinearities a, also satisfy 

3d>0, V* = 1,...,I, sup <n(r) <d. (6) 

rG[0,oo) J 



Proposition 2.1 Let T > 0. Assume (3), (6) and^i = 1,...,I, u® G L 2 (f2; [0, oo)). Then, there exists 
a nonnegative solution u = (ux, •••,«/) to i/ie system 

Mi = I,..., I, 

Ui G L 2 (Q T ), Jj ai(£0ui G L 2 (0, T; # 2 (ft)) , 



«, 



«; 



G L 2 (0, T; # 2 (^)) , fi» - 5iAui = itj on Q T , u { > (7) 

- A(/ ai(u)Ui) = u® on Q T , 



d n Ui = = d n (f Q ai(u)ui) on £ T . 

To prove Proposition 2.1, we will use the classical Leray-Schauder's approach, namely (see e.g. [9], 
Theorem 11.3) 

Lemma 2.2 (Leray-Schauder) Let (X, \\ ■ \\x) be a Banach space and T : X — )• X a continuous 
compact mapping. Suppose that 

3M > 0, Va G [0, 1], [ u G X, u = aTu ] => [ ||u|| * < Af ] . 

Then, there exists u £ X such that u = Tu. 

To define the mapping T, we will use the following lemma. 

Lemma 2.3 Let T > 0, w G L 2 (ft; [0, +oo)), ^ G L°°(Q T ), a, a £ (0, oo) suc/i too* < a < A < 
a < +oo. Then there exists a unique nonnegative solution w = w(A,wq) to 



w£L 2 (Q T ), Jj Aw £L 2 (0,T;H 2 (fi)), 

io — A ( J Aw J = uiq on Qt, <9 n ( J Aui J = on Et- 



(8) 



Moreover, if 

A n G L°°(Q T ), 0<a<A n <a<oo,A n ^A a.e., w% -»• w m L 2 (0), 

t/ien w(j4 n ,t(;g) converges strongly in L 2 (Qx) to w(A,wo). 

Proof of Lemma 2.3: Using convolution, we approximate A by a sequence of smooth functions 
(^4 n ) nG N G C°°(Qt) such that a < A™ < a and A n ^ A a.e.. Let also u>o be a regular approximation 
of u;o- There exists a classical regular nonnegative solution w n of (see e.g. [14], Theorem V.7.4, applied 
to the unknown A n w n ) 

d t w n - A(A n w n ) = 0on Q T , d n (A n w n ) = on E T) w n {0, •) = u#. (9) 



(10) 



Integrating (9) in time gives 

w n {t) - A ( f A n wA = < on Q T , d n ( f A n wA = on E T . 

We multiply by A n w n and use the following identity, valid for z n = A n w n : 

- f z n A f z n = f Vz n V f z n = f l-dt\V f z n \ 2 . (11) 

Jn Jo Jn Jo Jq 2 Jq 



We obtain the following estimate after integration in time 



f A n (w n ) 2 +fhvf A n w n \ 2 =[ w G l A n w n . (12) 

Jq t Jil 2 Jq Jq t 



In particular 

1/2 / r \ 1/2 



af K) 2 <aVf ( I «) 2> ) If (w n ) 2 ) ^a|K||^ (QT) <Zn/T|K 
Jq t \Jn J \Jqt J 



\L 2 (n)- 



(13) 



Now, up to a subsequence, w n converges weakly in L 2 (Qt) to some w. By the pointwise and 
uniformly bounded convergence of A n to A, for all ip G L 2 (Qt), tpA n converges strongly in L 2 (Qt) to 
•0^1 (using the dominated convergence theorem). Thus, Jq ipA n w n converges to J Q ipAw. In other 
words, z n = A n w n also converges weakly in L 2 (Qt) to z = Aw. 

By (10), A J" z n is bounded in L 2 (Qt); since j* z n is bounded in L 2 (Qt) as well, this implies that 
L z n is bounded in L 2 (0, T; H 2 (Q)). We now may pass to the weak limit in (10) to deduce that w is 
solution of (8). 

For the uniqueness, let w be the difference of two solutions of (8) (then w(0) = 0). We denote 
S(t) = f Q Aw. Formally, the idea is to multiply the equation w — AS = by S" = Aw. Then, after 
integration 

f Aw 2 =f S'AS = - f VS'VS = - f ld t \VS{t)\ 2 = - f l\VS{T)\ 2 <0. 

Jq t Jqt Jqt Jqt 2 J®. 2 

Whence w = since A > 0. 

Since we do not know whether VS" G L 2 (Qt), we have to justify this computation in an approximate 
way. For h G (0, T), let us denote 

V/iG(0,T), 5„(t) := g£±MzgW = i j/+ fe (^)( s )rf s . (14) 

Note that 

5 ft G L 2 (0, r - /»; H 2 (n)) , ||5 h - ^|| L 2 (Qt ^ -> as /» -> 0. (15) 

We have 

V« G [0, T - n), io(< + n) + w(t) - A [S(t) + 5(i + h)] = 0. 

We multiply by Sh(t) and integrate over Q to obtain 

/ [w(t + h) + w(t)]S h (t) = - f VS h (t)[VS(t + h) + V5(t)] = - f \ {\VS{t + h)\ 2 - \VS(t)\ 2 } . 
Jn Jn Jn n 



After integration on [0, T — h] and an easy change of variable, we have: 



[w(- + h) + w]S h 



To pass to the limit as h — > 0, we use 



(T-h,T)xU 



|V5| 2 + i [ |V5| 2 <i / |VS| 2 . 

h i(0,/i)x!i h JQ h 



(16) 



J \\7S\ 2 = [ -Sw= I - I w(t) I (Aw)(a)da dt < \\A\\ Laa{QT) h I w 2 dt 
JQh JQh Jn Jo l Jo J JQ h 



Now, letting h decrease to in (16) and using that 5/, — ^ Aw in 1? (see(15)), lead to Jq 2w Aw < 0, 
whence w = 0. 

Let us now prove the continuity result. Let us first notice that, for any solution of (8), we have the 
identity 



/ Aw 2 + / -|V / Aw\ 2 = \ w Aw 
Jq t Ju 2 Jo Jq t 



(17) 



This may be justified as we did above for the uniqueness (namely in the case wq = 0) by passing to 
the limit in the following identity where S(t) = L Aw, Sh(t) = [S(t + h) — S(t)]/h: 



i [w(- + h)+w]S h + VS h V[S(- + h) + S\=2 i 
Jo-r-h JO' 



I [w(- + h) + ui\S h + \ f 



(T-h,T)xQ 



\V Si- 



ft 



(o,h)xn 



\vs\- 



w S h , 



Qi 



w S h . 



(18) 
(19) 



And we pass to the limit as above as h — > to obtain (17) (at least a.e.T). 
Let w n = w(A n ,WQ). As in the beginning of this proof (see (17), (13)), the relation 



f A n (w n ) 2 + [ -|V [ A n w n \ 2 = [ w^A n w n . 
Jq t Jn 2 Jo Jq t 



(20) 



proves that w n is bounded in L 2 (Qt)- From (10), we deduce that L A n w n is bounded in L 2 (0, T; H 2 (Q)). 
A subsequence of (w n , A L A n w n ) converges weakly in L 2 {Qt) 2 to (w, A L Aw) and w is solution 
of the limit problem (8). By uniqueness, the full sequence converges. Since A n — > A a.e., y/A n w n 
converges also weakly in L 2 (Qt) to \jAw and, by the estimate (20), V J" A n w n converges weakly in 
L 2 (J7), the limit being necessarily VL Aw. In particular 



/ Aw 2 < liminf I A n (w n ) 2 , I |V / Aw\ 2 < liminf I \V I A n w 
Jq t n ^°° Jq t Jn Jo n ^°° Jn Jo 



™„.."|2 



(21) 



But, since lim^^co Jq WQA n w n = Jq wqAw, and since the identity (17) is true for w, it follows from 

(20), (17) that equality holds in the two inequalities (21). In particular, the norm of \/A n w n in L 2 (Qt) 
converges to the norm of \/Aw; this implies that the L 2 (QT)-weak convergence of \/A n w n to \/Aw is 
actually strong. Using again the pointwise convergence of A n , we deduce that w n converges strongly 
in L 2 (Qt) as well. 

□ 



Remark 2.4 As a consequence of (17), there is a constant C = C(a,a, H^oHi, 2 ^)) such that for any 
solution w of (8), 

\\w\\ L 2 (Qt) < VTC. (22) 

The next step is the definition of a compact continuous mapping T whose fixed points are solutions 
of (7). We introduce the Hilbert space 

X = UxKi^Xi, X l = {v€ L 2 (Q T ) : d t (J Si v) G L 2 (Q T )}, (23) 

where the Hilbert norm || • ||j is defined on X{ by 

Nl« : = IMIl2(q t ) + II^( j 5^)||!2 (Qt) , 

and where Jg = (I — SA)^ 1 is the resolvent of the Laplace operator on L 2 (ft) with homogeneous 
Neumann boundary conditions, that is 

[/ G L 2 (n), Z = J 5 f] <* [Z G H 2 (n), Z - 5AZ = f, d n Z = on dil]. 



Definition 2.5 We fix u° G L 2 (Q, [OjOo)) 7 . Let v = (vi,...,vj) G X and let u = (ui,...,uj) be the 
solution of (see [6], Proposition 9.24 an d Theorem 9.26): 

Vi = 1, ••■,/, Ui G L (0,T; H (O)), U{ — <5jAttj = Vi on Qt, d n Ui = on 'Ex- 
Next, we define 

T : X ->■ X by T(v) := u = (ui, ..., ui), 

where u% is the solution w of (8) with A = ai([u] + ),wo = u®; [u] + = ([u\] + , ..., [it/] + ) and [ui\ + is the 
positive part of Ui . 

Proposition 2.6 Assume (3), (6) and \/i = 1,...,J, it? G L 2 (Q; [0, oo)). Then the mapping T is 
continuous and compact from X into itself. 

Proof of Proposition 2.6: First, remark that for v G X, u = T(v) G X. Indeed, since Ui is 
solution of (8) with A = a,i([u] + ) and w$ = u®, we may write 

JSiUi = JsA I Au i + J k u °i = I AJ Si (Aui) + J Si u1 =>- d t (J Si Ui) = AJ Si (A Ui ) G L 2 (Q T ). 
Jo Jo 

Let v n be a bounded sequence in X. Up to a subsequence, me may assume that v" converges weakly 
to Vi in L 2 (Q T ). Then 

u? - 6iAu? = v? on Q T , d n v?} = on S T => d t v% = d t {J Si vf). 

Thus uf is bounded in L 2 (0,T;H 2 (n)) and d t u? = d t (J Sl v?) = J Si (d t vf) is bounded in L 2 (Q T ). As 
a consequence, -up is relatively compact in L 2 (Qt), and so is ["u™] + . Up to a subsequence again, we 
may assume that they converge strongly in L 2 (Qt) and a.e. in Qt- By continuity of a^, aj([tt n ] + ) 



converges a.e. and < d < aj([u n ] + ) < d < oo. By Lemma 2.3, u n := T(v n ) converges (up to a 
subsequence) strongly in L 2 (Qt)- Moreover 

< = A ( jf* a t ([u n }+) <) + u* => d t (J Si u?) = AJ 5i [( ai [u n }+) uf] . 

But the Yosida approximation AJs i is Lipschitz continuous on L 2 (Qt), and aj([ , u n ] + )u" converges in 
L 2 (Qt)- Therefore, ^(J^u™) converges also in L 2 (Qt)- Finally, this proves that u n converges in X 
(at least up to a subsequence) , whence the compactness of T. 
For the continuity of T, let v n — > v in X as n — > oo. If u n = (u™, ..., u 1 }) is the solution of 

Vi = 1, ...,/, u? - 5iAu™ = v? on Q T , d n u? = on E T , 

then £t" converges in L 2 (0,T; i7 2 (fi)) to the solution u% of 

ttj — <5jAuj = Vi on Qt, <9 n itj = on Sy. 

By definition, -u n = T(v n ) = (it", ■■■,u™) is the solution of 

f < GL 2 (Q T ), /;a,([i»]+)<a 2 (0,T;if 2 (!l)) 



fn ai([u n ]+)<) = n° on Q T , d n (/o Oi([« n ]+K) = on S T 



Using the compactness of 7" proven above, the sequence (u n ) ng N is relatively compact in X. Let 
■u°° = linip^oo u rLp be a limit point. Up to a subsequence, uf converges a.e. to ttj. By continuity of 
a,i, ai([u Up ] + ) — >• Ai := ai([u] + ) almost everywhere, and it is uniformly bounded from above and from 
below. According to Lemma 2.3, we can pass to the limit as n p — > +oo in (24). By the uniqueness 
result in Lemma 2.3 with A = Ai, we necessarily have u°° = T(v). The sequence (« n ) ne N lies in a 
compact set and has a unique possible limit point, so u n = T(v n ) — > T(v) and T is continuous on X. 

U 

Proof of Proposition 2.1: Let T G (0, oo) and a G [0,1]. Suppose that u G X is a solution of 
u = o~T(u) . By definition of T, we have 

( Vi = l,..,7,^£L 2 (Q T ), Ul >0, 

«,,/„* <*(«)«, G L 2 {0,T;H 2 {n)), 

Ui — 5iAiii = Ui on Qt, d n Ui = on Et, 
, m - A/ ai(u)ui = an,® on QT,d n (f ai(u)ui) = on E T . 

The initial conditions an? are uniformly bounded in L 2 (£l) for a G [0, 1]. Therefore, by the estimate 
(22), the function Uj remains bounded in L 2 (Qt), independently of a. We also have dt{Js t Ui) = 
AJs t (ai(u)ui), so ii is bounded in X independently of a. Using Proposition 2.6 and Leray-Schauder's 
Lemma 2.2, we can conclude that T has a fixed point, which is a nonnegative solution of (7) (the 
nonnegativity of u, is a consequence of u% > and of the maximum principle property of (7 — (5jA) _1 
with homogeneous Neumann boundary conditions, see e.g. [6], Proposition 9.30). 

3 L°°-estimate of u in Proposition 2.1 

A main estimate in the proof of Theorem 1.1 is given in the next proposition. 



Proposition 3.1 Assume u° £ L°°(0, [0, +oo)) J and (3), (6) as in Proposition 2.1. Let us define 

V/c>0, G(/c)=max{ sup a^r)}. (26) 

* rG[0,fe] 7 

Then, for any solution u,u of Proposition 2.1, we have 

max \ 5i\\ui\\L°°(Q T ) + II / Oi(«)wi||z,°°(Q T ) f - M o + MiTG(k ), (27) 

where Mq, M\ and k$ depend only on u°,5 := minj <5j, o" := maxj Si. 
The proof of Proposition 3.1 uses the following classical lemma. 
Lemma 3.2 Let f £ L°°(r2) and Ze£ w satisfy 

w £ H 2 {yi), w > 0, —Aw < f on £1, d n w = on t?$7. 
T/ien i/iere exists C = C(£l) suc/i that 

IMU«(n) < c (ll/IU~(n) + / w ) • ( 28 ) 

Proof: First, we rewrite the equation as if — Aw; < / + to. Let us fix p £ (N/2, oo). Using u; > 0, the 
comparison principle and elliptic regularity theory, we know (see e.g. [9], Theorem 8.15) the existence 
of C = C(Vt,p) such that 

IMIioc < C(||/ + «;|| L p) < C(\\f\\ LP + HwIIlp) , 

< C (\\f\\ LP + IMlfc 1)/p (/^) 1/p ) , 

^ C I ll/lli> + elMli 00 + c ( £ ) / w ) (Young's inequality) 

and we conclude choosing e small enough. 

□ 

Remark 3.3 Obviously, the conclusion of Lemma 3.2 would be the same when assuming only / £ 

LP{n),p> N/2. 

Proof of Proposition 3.1: We rewrite the equations in Ui,Ui of Proposition 2.1 as 

Ui-AySiUi+j ai{u)uA = u® , Ui-Awi = u®, Wi = 5^ + / a* (&)«». 



(29) 



We sum up the equations (29), denoting U = J2i Ui,W = J2i w i '■ 

U-AW = U°:=J2 u i- ( 30 ) 

i 

Next, we apply Lemma 3.2 with w = W(t),a.e.t,f = U° (note that -AW(t) < U°). It gives 

a.e.t, \\W(t)\\ Laa(n) < C (\\U°\\ L oo {U) + J W(t)j . (31) 

9 



By nonnegativity of Uj, ai(u)ui, we also have (see the definitions of W, Wi): Vi = 1, ...,/, a.e.t G [0, T\: 

5i||^(t)|| L oo (Q) ,|| / ai(u)Ui\\ L oo(ty < \\W(t)\\ L oo(ay 
Jo 

Then, to end the proof of Proposition 3.1, it is sufficient to prove the following lemma. 

□ 

Lemma 3.4 

a.e.t G [0, T], [ W(t) < C + dTGfa), 

where Co,C\,ko depend only on u°,8,5 and G is defined in (26). 

Proof of Lemma 3.4: By integrating the equations on m and «j in Proposition 2.1, we get: 

Vt > 0, f Ui(t) = I Ui(t) = [ u°. (32) 

Jn Jn Jn 

Recall that Ui,Wi G L 2 (0,T;H 2 (£l)),ai(u)ui G L 2 (Qt). We also have d t Ui = AJs^a^ujUi) G L 2 (Qt) 
From (29), we may write with dtWi = SidtUi + ai{u)ui G L 2 {Qt) 

d t Wi - 5iA(d t Wi) = ai{u)ui. (33) 

Differentiating d n Wi = with respect to t on dtt leads formally to d n (dtWi) = 0. Let us check that 
dtWi = 6{t) where 9(t) is the unique solution of 

9eL 2 (0,T;H 2 (n)), a.e.t G[0,T], 9(t) - 6iA6(t) = (oi(fi)«i)(t), d n 6{t) = on dtt. (34) 
Using also ai(u)v,i > 0, it will then follow that 

d t Wi>0, d t Wi £ L 2 (0,T;H 2 (Q)), \\d t Wi\\ L 2(Q T ) < \\ai(u)ui\\ L 2(Q T y (35) 

By integration in time of (34), and with 0(£) = J 6(s)ds, we have 

8(i) - 5 l AQ(t) = I (oi(u)ui)(t), d n Q(t) = on dtt. 



Comparing with Wi — 5iAwi = 5^ + J ai(u)ui, d n Wi = implies by uniqueness that: 
0(i) = Wi + {I — SiA)^ 1 ^, whence Q'(t) = 6 = w\ after differentiating in t. 

We denote 

V = ^2 SSi, B =} j ai(u)ui. 



Recall also that 



U = ^2ui,W = ^2wi,Wi = SiUi+ ai(u)ui. 



Summing the / equations in Ui,Ui as in (30), we have 

S~ 1 V-AW<U-AW = U°. (36) 

10 



We multiply this equation by dtW = J2i dtWi = &tV + B > (see (35)) and we get 



-?-i 



1 



5 / V(d t V + B) + - / d t \VW\ z < / U u (d t V + B). 
Jn * Jn Jn 

We integrate in time to obtain (we denote V° := V(0) = W(0)) 



V\T)+ / 2BV + 5 / \VW{T)Y < / (V") 2 + 5\VV U \ 2 + 25U U {V(T) -V iJ ) + / 25U"B. (37) 
a Jq t Jn Jn Jq t 



Since we have by definition 

5U° = 5U° - 5AV° >V°- 5AV°, 

we have 

/ {V ) 2 + 5\VV Q \ 2 - 25U°V° < - [ (V ) 2 + 5\VV°\ 2 < 
Jn Jn 

So that (37) becomes, 

/ V 2 (T) + I 2BV + 5 I \VW(T)\ 2 < 25 I U°V(T) + I 25U°B. 
Jn Jq t Jn Jn Jq t 

We have in particular, with ||f7 ||oo = ||£^ ||l°°(T2)> and by using (32): 

I BV<~5\\U Q \\ 00 \ 

Jqt 

Thus, we have for any k > 



s / B < 5\\U°\ 

lQ T n{v>k} 



n JQ 



l V°+ _ B+ _ B 

n JQ T n{v<k} JQ T n{v>k} 



Note that, {V < k} C C\i{ui < k5 }. Thanks to the L 1 estimate (32), we have 



/ 

JQ- 



B 



[ _ y;ai(fi)«i<r /[/ 

iQ T n{y<fci , un 



G^ 1 ), 



? T n{v<fc} 7Q T n{y<fe} 

where G is defined in (26). Finally choosing k = 2<5||£/ || oo in (40), we obtain 



/ - 

JQ T n{V 
Adding the two last inequalities gives 



B < ( 2 / V° + T 

T n{v>k} V Jn 



U [ 



irx\\TT0\ 



G(k ) , k = 26r l (5\\U 



f 5<C + C 1 rG(fe ) 



where Ci depends only on u°, 5,5. 



To end the proof of Lemma 3.4, we use that W{t) = ^ 5iUi(t) + j* B(s)ds so that 



Vie [0,T], / W(t) < / U°+ / 5. 

Jf7 Jf7 ./Q T 



(38) 



(39) 



(40) 



(41) 
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D 
From the L°°-estimate of Proposition 3.1, we may now deduce that the problem (4) in Theorem 1.1 
has at least a weak solution under the only assumption of continuity of the a^'s. 

Corollary 3.5 Assume (3) (only) andMi = l,...,I,u® G L°°(Q; [0, oo)). Then, there exists a nonneg- 
ative solution u = (in, ...,ui) to the system 

' VT>0, Vt= 1, ...,/, 

Ui,Oi{u)ui G L 2 (Qt), /oOi(«K G £ 2 (0,T;tf 2 (ft)) , 

«i € £°°(Qt) n L 2 (0, T; F 2 (0)) , fij - SiAiii = Ui on Q T , (42) 

m - A(/ ai(u)ui) = u® on Q T , 

d n Ui = = «9 n (/ ai(u)ui) on T, T . 

Proof: Here, aj is assumed to satisfy only (3) (and not (6)). Let T > 0. We introduce M 2 := 
cp 1 [Mo + MiTG(ko)] where the function G is defined in (26) of Proposition 3.1 and M ,Mi,k are 
defined in (27) of the same proposition. We define 

Vr G [0,M 2 ] 7 , ai(r) := ai {r),Vr G [0,oo) 7 \ [0, M 2 ] 7 , a,(r) = min{ a 4 (r), G(M 2 )}. 

Then, aj is continuous on [0, oo) 7 and 

< d < a~i < G(M 2 ) < oo, 5j < a-/. 

Therefore, we may apply Proposition 2.1 with a, replaced by 5j. By Proposition 3.1, the corresponding 
u satisfies 

Vt = 1,..,/, ||«i||Loo (QT) < r 1 [M + M l TG(k )) , 

where G is defined as in (26) with a* replaced by Sj. But G(ko) < G(ko), so that 

Vi = 1, ..., /, < Uj < M 2 , ai(tt) = aj(ix). 

Therefore, the solution obtained with the data a, is also solution with the data ctj. 

This provides a solution of (42) in Corollary 3.5 with the estimate (27), but only on [0, T] and it 
may depend on T. To construct a global solution on (0, oo), we may argue as follows: let T p be an 
increasing sequence of times with limp^. +00 T p = +oo. Let u p be a solution of our problem on the 
interval [0, T p ] given by the above proof. For k G N, we denote by X k the space X as defined in (23) 
with T replaced by T^ and we denote by T k : X k — > X k the operator T with T = Tj.. For p > k, we 
denote u P)k := u? 0T , so that T k (u p,k ) = u p ' k . We will prove that 

\/k G N, (u p,k ) p >k is relatively compact in X fc . (43) 

Thus, using a diagonal process, we obtain a sequence p m — > oo as m — > oo and some limit u defined 
on (0, oo) so that, for all k G N, u Pm,k converges to M[o,r fc ] m X k as m — > oo. Then, 7fc(«ro,T fe ]) = u [o,T k ] 
and it is a global solution of (42). 
Let k be fixed in N and let us prove (43). By the L°°-estimate (27) in Proposition 3.1, 

Vp > k, \\u p \\l°°(q Tl) < Um + M.nGiko)}. (44) 

Thus, ai(u p ) is uniformly bounded on Qx k - This implies by (22) that u p is bounded in L 2 {QT k Y and 
so is dtu p since by (35) 

h\\dtu p \\ L 2 {QTk) < 2||a 4 (-u p )nf|| L 2 ( Q Tfc) < C{k). 

Thus, u p,k is bounded in X k and, by compactness of T k , it is relatively compact in X k , whence (43). 

□ 
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4 Proof of existence in Theorem 1.1 

Existence of a weak solution to (4) is already proved in Corollary 3.5. It only remains to prove that 
this solution is actually as regular as stated in Theorem 1.1. This will mainly be a consequence of the 
L°°-estimate on u proved in the previous section, namely 

Vz = l,...,I, ||«i||ioo(Q T ) < Cb + CiT, \\ai(u)\\ Lx{ Q T) <C(T), 

where Co, C\ depend only on the data and C(T) = G(C$ + C\T). 
We begin by the following simple estimates. 

Proposition 4.1 Let Wi = 5iUi + f Q ai(u)ui where u,u is solution of (4.2) in Corollary 3.5. Assume 
u° GL°°(ft,[0,oo)) 7 . Then, 

VT > 0, Vw t G L°°(Q T ) N , Wi, d tWi G L°°(Q T ), d tWi > 0. (45) 

Proof: The fact that Wi G L°°{Qt) is a consequence of (31) and Lemma 3.4. We recall the two 
equations (see (29), (33)): 

vn - Awi = u°, d t Wi - 5iA(d t Wi) = ai(u)ui. 

Since Wi, Awi G L°°(Qt) and d n Wi = on Sy, we deduce that Vwi G L°°(Qt) n (at least). We have 
already seen that dtWi > comes directly from the second equation and the nonnegativity of ai{u)ui. 
Now we rewrite this equation as 

(d tWi - C{T)ui) - 5iA{d tWi - C{T)ui) = ( ai (u) - C{T)) Ui < 0. 

Together with d n (dtWi — C{T)ui) = on St, this implies 

d t Wi - C{T) Ui < 0, so that < d tWi < C(T)[C + dT]. 

□ 

We will now prove that Ui(t,x) := f Q [ai(u)ui](s, x)ds is in C a (Qr) so that, since U{ = Wi — Ui, it 
will follow that Ui is not only bounded, but Holder-continuous (at least). 

To prove it, we rely on the C a -regularity theory of Krylov-Safonov for the solutions of nondivergence 
parabolic equations with bounded coefficients. We actually use them in the rather particular case of 
the operator — AA where A is bounded from above and from below. We may state the result we need 
as follows: 

Lemma 4.2 Let A G C{Qt),9 G L oc (Q T ),a,a G (0,oo) with 0<a<A<a<oo. Let w G 
C 2 ' 1 (Qr)nC 1 ' 1 (Qr) solution of 

d t w - AAw = ginQ T . . 

d n w = on E T , MO) = 0. l ' 



Then, there exists a G (0, 1), C > such that 



IT 
where a,C depend only on a, a, T, \\g\\ L <x>(Q T \,Q. 



w\\P < C (47) 
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Remark 4.3 Note that the L°° estimate of w is easy by a comparison argument (valid here thanks to 
the a priori regularity of w and of A ) : we remark that the function W(t, x) := t sup g is a supersolution 
of the problem (46), so that W > w. Doing the same from below, we obtain 

IMIl°°(q t ) < r||^|| L oo(Q T ). (48) 

Next, we may use the Krylov-Safonov result: the global estimate with homogeneous Neumann bound- 
ary conditions as stated above may, for instance, be found in [10], Lemma 2.2 (in a quite more general 
setting). We more generally refer to [13], [10], [20] for this kind of results. 

We apply this result to prove the regularity of Ui = f Q ai{u)ui. 

Proposition 4.4 Let T > and u° G L°°(0, [0, oo)) 7 . There exists a G (0, 1),C > such that 

ii/7.||( Q ) _l iif,.i|( a ) < r 

\\ u i\\t ' II u iIIt — 

Proof: Let u,u be the solution of (42) in Corollary 3.5. Recall that < d < ai{u) < C(T). Since 
Lemma 4.2 a priori applies to regular solutions only, we will use a convenient approximation of u. For 
this, let A n be a smooth approximation of a,i(u) such that 

0<d<A n <C{T), A n ^ ai {u)a.e. 

Let also v n be a smooth approximation of u® such that 

< v n < ||«? \\ L oo {a) , v n -> u° in L 2 (n). 

Let uf be the solution of 

d t uf - A(A n uf) = 0, d n u™ = on S T) <(0) = t; 71 . 

Then, after integration in time, we see that Uf 1 = J A n u™ satisfies 

8tU? - A n AU? = A n v n , d n U? = 0on E T , U? (0) = 0. (49) 

By Lemma 4.2, there exists a,C independent of n such that ||t/"||j* < C. By Lemma 2.3, uf 
converges to Ui in L 2 (Qt) which implies that U- 1 also converges to Ui in L 2 (Qt)- Whence the estimate 
of Proposition 4.4 on Ui. The estimate on U{ = W{ — U- h follows by combining with Proposition 4.1 
which says that Wi is even Lipschitz continuous. 

□ 

Now that we know that the coefficient ai(u) is not only bounded but also continuous, we may 
continue improving the regularity of u. 

Proposition 4.5 Assume u° € L°°(0, [0,oo)) 7 . Then, 

Vp G [1, oo), VT > 0, V? = 1, ...,/, Ui , d t Ui, AU G L P (Q T ). 
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Proof: We may formally write 

dtUi - ai {u)AUi = ai(u)u1, d n Ui = 0, tfi(0) = 0. (50) 

Here ai{u) is continuous on Qt so that, aj(n) being given, this equation has a unique solution : let us 
call it V{. We set V{ := dtVi/di(u). Then 

t 



AVi = u°, Vi= di(u)vi, d n { / ai(u)vi) = 
Jo Jo 



Thus, Vi coincides with our u% (and Vi coincides with our Ui) thanks to the uniqueness result of Lemma 
2.3. 

Moreover, L p -maximal regularity holds for the equation (50) since cii(u) is continuous (see for in- 
stance [14], Theorem 9.1 or [21], Theorem 2.5.2) so that, as ai(n)n? G L°°(Qt) C L p (Q t ), we have 

VpG (l,oo), ||d t [/i|| LP( Q T) ,||A[/i|| LP( Q T) < C, 

where C depends on p, ||ai(w)tt°||^oo(g T ) and on the modulus of continuity of the function ai{u). 

Next, from < Ui < d~ ai{u)ui = dT \dtUi\, we deduce that Ui G L p (Qt) as well. And p = 1 is also 
included since Qt is bounded. 

□ 

With Proposition 4.5, the first part of the existence result in Theorem 1.1 is now complete. We will 
now assume that a% is locally Lipschitz continuous. 

Proposition 4.6 Besides (3), assume a% is locally Lipschitz continuous for alii = 1,...,I. Assume 
also u° e L°°(n, [0,oo)) 7 . Then 

Vi = l, ...,/, VT>0, m eL°°(Q T ), Vp<oo,Vr G (0,T), d t u h A(ai(fi)«i) G L p ((t,T) x 0), 

and 

djiij - A(ai(u)ui) = on Q T , d n (a,i(u)Ui) = on E T , 

is satisfied pointwise. 

Proof: The equation in Uj may also be written (at least formally to start): 

d t (ai(u)ui) - ai(u)A(ai(u)Ui) = UiDai(u) ■ d t u. (51) 

We know that 5idtUi + Oi(u)ui G L°°(Qt) (see Proposition 4.1), and ai(u)Ui G L p (Qt), for all p < oo, 
so that dtii-i G L p (Qt) for all p < oo. The right hand side of this equation F := UiDai{u) ■ dtu is 
therefore in L p (Qt) for all p < oo since also Da,i(u) G L°°{Qt) (as a-i is locally Lipschitz and u is 
bounded). 

Since ai(u) is continuous on Qt, we know (see again e.g. [14], Theorem 9.1 or [21], Theorem 2.5.2), 
there exists a (unique) solution 9 to 



Vp<oo, 9e C(0,T;LP(O)),VrG (0,T), d t 6, AO G L p ((t,T) xfi) 
<9 t # - ai(ti)A# = F, d n = on S T , 0(0) = ai(n°)n° 



o^„o (52) 



and we have 

I|0||l°°(Qt) + \\dt®\\Lp((T,T)xn) + \\A9\\LP(( T ,T)xn) < C[||-^IIlp(q t ) + IK IU°°(n)]j ( 53 ) 
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where C depends on r, T,p, Q and of the modulus of continuity of a,i(u). 
If we knew that 9 = ai(u)ui, then the proof of Proposition 4.6 would be complete using moreover: 

d t Ui = ai(u)~ l [d t {ai(u)ui) - UiDai(u) • d t u] G L p ((t,T) x O). 

To prove it, we recall (see the proof of Lemma 2.3) that m is the limit of the approximate solutions 
u n of 

d t u n - A(A n u n ) = 0, d n u n = on £ T , u n (0) = uf, 

where A n is smooth and converges pointwise to a,i(u) with < minaj(u) < A n < maxaj(u) < +oo. 
Moreover, u n is bounded in LP(QtY for all p < oo by the analysis in Proposition 4.5. Here, we choose 
such an approximation A n which moreover satisfies 

A n -> oi(fi) m L°°(Q T ), d t A n -> d t a;(u) = Doi(it) • d*u in L P (Q T ) Vp < oo. 

Then, we apply the estimates (53) to A n u n which satisfies 

d t {A n u n ) - A n A{A n u n ) = u n d t A n , d n A n = on S T , A n u n {0) = Oi(t2°)«?, 

and they are preserved at the limit. Whence = a,i(u)ui by uniqueness in (52). 

□ 

Proof of the Existence in Theorem 1.1: It is a consequence of Corollary 3.5 and of Propositions 
4.4, 4.5, 4.6. 

□ 

Remark 4.7 Note that, not only we proved existence of a solution with the announced regularity, 
but we even proved that any weak solution as in Corollary 3.5 has actually the announced regularity. 
This will be useful in the proof of uniqueness 

Remark 4.8 The assumption that the ctj are bounded from below is essential in our proof of existence, 
first for the L 2 -estimate, next to apply the Krylov-Safonov regularity theory. In the case when the 
cii degenerate (oj > 0), the L 2 a priori estimate is to be replaced by yja,i{u)ui G L 2 (Qt)- However, 
we loose the L 2 -compactness of the approximate solutions and also most regularity properties of the 
solution as well. It could however be interesting to study the possibility of existence of weak solutions. 

Remark 4.9 The above analysis relies on the use of the C Q -Krylov-Safonov estimates. However, it 
is interesting to notice that one can prove directly, by an elementary estimate, that u G L°°(Qt), 
without using these estimates in the (rather general situation) where, besides (3), a, satisfies 

Vi = 1, ...,/, aj is locally Lipschitz continuous, Vj = 1, ..., I, 9a -Oj > 0. (54) 

Once the L°°-estimate is proved on m, the full regularity follows by the same arguments as in Propo- 
sition 4.6. 
We indicate below (at least formally) the computations which leads to u G L°°(Qt)- 

Proof of u G L°°(Qt) under assumption (54): 

We write a^ for a,i(u) and aij = d^-di. We multiply the equation dtUi — A(ajUj) = by p(aiUi) p ~ 1 
and we integrate over fi: 



f af-y + f p(p- l)(a lUi y- 2 \V(a lUl )\ 2 = (p - 1) V / a\- 2 u\a l3 d t u r (55) 



Tt> 
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We proved in Proposition 4.1 that dtUj + cljUj < C(T) < oo. This implies dtUj < C(T). Plugging this 
into (55), using a%j > 0, aij bounded and en > d leads with some Ct independent of p to: 



d 
~dt 



I a?" 1 "? + f P(P~ l)(am) p - 2 \V(a lUi )\ 2 < C T (p - 1) V / a?" 1 ^ . 
Summing over i and using Gronwall's lemma on the term £^ J*„ a^~ n^, we then have 

E/arvw<e 7T ^^E/«rV(o)- 

Using the lower and upper bounds on aj, we have with ^4, C\ both independent of p: 

£ / «(t) < Ae c ^(l + £ HcviiCOJIIooy. 

This implies 



a,; it. 



)(t)|| P <^ 1/p e C T(l + ^|| a ^( )|| c 



whence the L°°-estimate on aiUi by letting p — >■ oo, and then on Uj itself by using the lower bound on 

□ 

5 Proof of uniqueness in Theorem 1.1 

Actually, we will prove the more general following result: 

Proposition 5.1 Let u E L°°(Q, [0, oo)) . Assume that for all i = 1, ...,/, a% satisfies (3) and is 
locally Lipschitz continuous. Then there exists a unique solution to the system (42) in Corollary 3.5. 

Proof: By Remark 4.7, we already know that any solution of (42) satisfies the regularity stated in 
Proposition 4.6 and Theorem 1.1. Let u,v be two such solutions. We denote a, = ai(u),bi = a«(u). 
By difference, 

d t (ui - Vi) - A [a,i(ui - Vi) + Vi(ai - &,)] = 0. 

We set 

U{ = Ui — Vi, Ui = Ui — Vi,U = u — v, Ai = J Da,i(tu + (1 — t)v)dt, 

Jo 

so that ai — bi = Ai ■ (u — v) = ^ ■ AijUj. Note that ||-Ai||l°o < oo. Then 



d t Ui - A 



aiUi + ViAi • U 



0, d v (a i U i + v i A i -U) = 0. (56) 



Lemma 5.2 Let F 6 C^ c (Qt) 1 ■ There exists a solution to the dual problem 



Vi = 1, ...,/, fi,d t fi,Aipie L 2 (Q T ), 

dt<fi + ajA^j + J^-Bj • A(/j) = Fj on Q T , (57) 

93 = (931, ..., ipi), d u ipi = on S T , <fi{T) = 0, 



w/tere 2?j = (Bn, ..., Bu), B i:j = VjAji. 
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Assuming this lemma, we multiply each equation (56) by (fi and we obtain after integration on Qt (the 
integrations by parts are allowed, thanks to the regularity of u, v, u, v, (fi and the boundary conditions; 
we also use J Qt UiJ &i (Bi ■ Aip) = J Qt U^Bi ■ Aip): 

0= / U i [d m + a l Ai Pi ] + A(p i v i A i -U= / U t Fi - UiBi ■ Aip + A^ Vi Ai ■ U. 

JQt JQt 

Summing these / identities gives Y2i Jn UiFi = which implies U = by arbitrarity of the Fi, whence 

uniqueness. 

□ 

Proof of Lemma 5.2: To solve the dual problem (actually interesting for itself), we may start with 
Oj replaced by regular approximations A™ converging in the usual way to at (which means a.e. and 
uniformly bounded from above and from below), and we first solve 

d t 0? + A(A n 0?) + A J*. (A • 6 n ) = AF U d n (A?9?) = 0, 9i(T) = 0. 

This is possible since 9 G L 2 (Qt) — >■ (AJs t (Bi • 9)) l<i<1 G L 2 (Qt) is a Lipschitz perturbation (recall 
that Bi G L°° and A J5. is the Yosida approximation of the operator —A with homogeneous Neumann 
boundary conditions). Note that J„ Of it) = 0. Next, we solve 

A99™ = of in n, d n (<p?) = m an, [ $ = o, 

Jn 
so that, "by applying A -1 " to the equation in Of, we obtain 

d t <pf + AfAipf + J 5i {Bi ■ Acp n ) = F, d n (<pf) = on S T , <p?(T) = 0. (58) 

Next, multiplying by Aipf gives 

[ -ld t \V<p?\ 2 + A?(A<p?) 2 + A<p?J Si {B i -A<p n )= [ FiAtf< [ e(Atf) 2 + C € F 2 . 
Jn * Jn Jn 

We choose e := d/2 and we deduce 

/ -k|V^| 2 + =(A^™) 2 < C / Ff+ I VZVpf <C f Ff+ f e\VZ\ 2 + C e \Vipf\ 2 , (59) 
Jn z * Jn Jn Jn Jn 

where Z - biAZ = Bi ■ Aip n , d n Z = 0. Multipling this by Z gives 

I Z 2 + 5i\VZ\ 2 = I ZBi- At/ 1 < ||5i|| £ oo I eZ 2 + C e \A<p n \ 2 =► / |VZ| 2 < C / 
Jn Jn in Jn Jn 

Summing the equations in (59) and choosing adequately e leads to (with a different C) 

-dt f £ |v^| 2 + = £ / ( A ^) 2 < c f / £[i? + |Vtf 

Jn i z . Jn Jn , 

Integrating the Gronwall estimate in £\ |V(/?™| 2 and plugging back the terms in A<pf yield 

sup /^|V^| 2 + =/ \A^\ 2 <c[ \F\ 2 . 

0<t<TJn ■ z JQt JQt 



\Ay 
'n 



.ni2i 



By going back to (58), we also obtain that dt<pf is bounded in L 2 (Qt)- Now, we can pass to the limit 
as n — > 00, weakly in L 2 (Qt) in each term of (58), to prove the existence result of Lemma 5.2. 
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□ 

Remark 5.3 We do not know whether uniqueness holds or not without assuming Lipschitz continuity 
of the a,i . The above proof indicates that uniqueness is essentially equivalent to solving the " dual" 
problem (57). The fact that Bi £ L°°{Qt) (which equivalent to the Lipschitz continuity of Oj) is 
strongly used in the estimates to solve (57). It is not clear how to weaken it. 

6 A constructive approximation procedure 

In this Section, we just mention without proof a complementary approach for the proof of existence 
and regularity of solutions to our system. A interesting point is that it provides a constructive approach 
and may be used to provide a numerical approximation schemes of the solutions (and have actually 
been used in [15]). 

The approximation procedure follows the ideas of [15] to make a time semi-discretization with an 
explicit treatment of u and an implicit treatment of Ui in ai{u)ui. We fix T > in the following. Let 
r > is given. We introduce the following approximate system 

in_l_1 _ 

,n+l 



S-^ - Aa i [(fi ft K +1 ] = 0, in ft, 

-diAuf + uf = <, in U, (60) 

d n u™ = d n u™ +1 = 0, on dn. 

Then, using some adequate discretized version of the approaches in Section 3, it may be proved that: 

- there exists a unique solution to the system (60), 

- the discrete version of the L 1 -estimate (41) is valid for B n = ^ ai(u n )u™ , 

- an L°°-estimate may be deduced for u n , 

- the direct L°° estimate may be reproduced in this direct situation under the assumption (54). Once 
all these estimates are obtained at the discrete level, one can get Corollary 3.5. 

A common issue to both procedures is the need, first to truncate the coefficients (for steps 2 and 3), 
and then notice a posteriori that this truncation was useless. 

Remark 6.1 We have used here a situation where the i-th relaxation parameter depends on the 
species (in words Si depends on i). It is quite natural on a modeling point of view to consider: 

• situations where the relaxation parameter might depend on the species that are observed and 
on the species that are observing (in words, ai = ai (J5 n (ui), ■ ■ ■ , J&uWi)))- 

• several sensitivities for a single species might have to be considered, short and long distances 
could matter in the diffusion coefficient (in words two or more different Jg(uj) could be taken 
into consideration in <Zj). 

Our results could be easily extended to these more general situations through rather slight modifi- 
cations. 
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